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Using 226 million BB events recorded on the T(4S) resonance with the BABAR detector at the 
SLAC e + e~ storage rings PEP-II, we reconstruct B~ — > D*°e~u e decays using the decay chain 
D*° — > D°ir° and D° — > K~tt + . From the dependence of their differential rate on w, the product of 



4 



the four- velocities of B and D*° , and using the description of the form factor F(w) by Caprini et 
al., we obtain the preliminary results p Al = 1.15 ± 0.06 ± 0.08, F(l)-\V cb \ = (36.3±0.6± 1.4) • 10~ 3 , 
and B (B~ -> D*°e"77 e ) = (5.71 ± 0.08 ± 0.41)%. The first errors are statistical and the second 
ones are systematic. 



Submitted to the 2007 Europhysics Conference on High Energy Physics, 
Manchester, England. 
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The exclusive B-meson decay modes with the highest 
rates are the two scmilcptonic modes B° — > D* + e~~P e 
and B~ —> Z?*°e~77 e . Whereas the first has been mea- 
sured by many experiments [l[ to determine its rate T, 
its differential rate dT/dw, and the CKM matrix clement 
iVciJ, the second has only been measured by two groups 
|2j, lla] with lower statistics. In the £>° mode, the ob- 
served differential decay rate is well described by heavy- 
quark effective QCD (HQET) using form factors with the 
slope parameter p 2 . However, the B° experiments do not 
agree well in their p 2 results. Using the isospin symme- 
try dT(B~ -► D*°e-V e ) = dT(B° -► D* + e~V e ), a preci- 
sion measurement for the B~ mode can help to improve 
knowledge of p 2 and consequently of T and \V c b\- 

The aim of our analysis [H is the determination of 
the differential decay fraction dB/dw(B~ — > Z?*°e _ F e ), 



where B is related to the decay rate V by the known 
lifetime t(B~) and w is the invariant product of the 
four- velocities of B~ and D*°. The neutrino in the 
B~ — ► D* a e~V e decay is not reconstructed. Therefore, 
the w value of each reconstructed event cannot be ob- 
tained, only an approximation w which will be defined 
below. Instead of unfolding, the parametrized dB/dw 
expectation together with the w resolution from Monte 
Carlo simulation (MC) is fitted to the observed dB/dw 
distribution. Our fit, as in other recent B° — > D*iv 
analyses, uses the parametrization of Caprini et al. [B| 
and determines the two parameters F(l) ■ \V c b\ and p 2 . 
The third result, B(B~ — ► D* e~F e ), is obtained by in- 
tegrating dB/w. The parametrization is defined by the 
following expressions: 



dr (B~ -> D*°e~V e ) Gj 



dw 



F(w) 2 = \h Al (w)\ 2 



487T 3 



(rriB — mrj*) 2 m 3 D , \/ w 2 — 1 (w + l) 2 



4w m 2 B — 2wiriB'mD* + n^n* 
+ 1 (m B - m D * ) 2 



F{w) 2 \V cb \ 2 , 



Aw m 2 B — 2wrriB'mD* + Tni- 



D* 



H (w) 



w — 1 , 
! + — {!- R 2 {w)) 



+- 1 


(m B - 


m D * 


-i 2 
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H±(w) 




- 15) z 2 


- (231p 2 Al - 



E 

i=0, + ,- 



\Hi(w) 



1 — 2wr 



VwTT- V2 



= i?i(l) - 0.12(w - 1) + 0.05(w - l) 2 , R 2 (w) = R 2 (l) + O.U(w - 1) - 0.06(w - l) 2 . 



r 



Note that p 2 = p\ in this notation. The parameters 
i?i(l) and R 2 (X) are not determined in this analysis; we 
use the BaBAR results from the B° — > D*lv decay as 
input, see @ and Table HI 



For our analysis, we use 205 fb 1 of e + e annihila- 
tion data recorded at y/s ~ m(T(4S)) with the BaBAR 
detector Q at the SLAC storage rings PEP-II ||. In 
addition to these on-peak data, we also use 16 fb -1 of 
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off-peak data collected below the T(4S) resonance. We 
select B~ — ► D*°e~17 e candidates Q by pairing electrons 
with p* > 1.2 GcV/c in the e + e~ rest frame with D*° can- 
didates. Since the precision of our results is not limited 
by statistics, we restrict the analysis to the sequential de- 
cay modes D*° — > £>°7r° and D° — > A"~7r + , which have 
the smallest combinatorial background and the best res- 
olution in Am = m(K~ir + ir°) — m(K~ir + ). 

Charged particles are selected if they have at least 
10 hits in the BABAR drift chamber, transverse momen- 
tum pt > O.lGcV/c, and a polar angle between 23.5° 
and 145.5°. Electrons (kaons) are selected with tight 
(loose) BABAR particle identification criteria. Neutral 
pions are reconstructed from two photons with energy 
above 30MeV and a photon-compatible lateral shower 
shape in the BABAR calorimeter. The invariant mass must 
be 115 < ?ti 77 < 150McV/c 2 , and the photon pair is then 
constrained to a common vertex and to m 77 = m(7r ). 
The decay candidates have to fulfill the following fur- 
ther requirements: The D*°-D° mass difference must be 
135 < Am < 153MeV/c 2 and the D°-candidate mass 
1.8496 < m(K--K+) < 1.8796 GeV/c 2 . To reject non-B- 
decay candidates, the normalized Fox- Wolfram moment 
i?2 [10| of the event has to be smaller than 0.45. To reject 
candidates with a D*° from one B meson and an electron 
from the other B in the event, the angle between the D*° 
and the e~ has to be larger than 90°. 

Since there are many low-energy background photons, 
the selection criteria result in many events with two or 
more D*°e candidates, on average 1.75 per event. All 
D*°e candidates in the same eKir combination are col- 
lected into one candidate group; on average there are 
1.015 candidate groups per event. Only one candidate 
group per event is kept, in case of multiple groups the one 
with the smallest deviation \m(Kiv) — m(D°)\. All can- 
didates in one group are kept in the analysis because the 
simulation of low-energy photons is not perfect. This pro- 
cedure ensures that correctly reconstructed candidates 
are selected with the same probability in data and MC. 

The set of surviving candidates is binned in three di- 
mensions according to their values of Am, cos#b Y , an d 
to. The first two variables are used for the separation be- 
tween signal and background, the third is used for the w 
dependence of the signal. Am is defined above, and 9 BY 
is the angle between the directions of the B meson and 
the Y = D + e system in the e + e~ rest frame under 
the hypothesis that the B decays into only D*°, e, and 
neutrino. It is defined by the four-vector relation 

pi = = ml + m\ - 2(E B E Y - \p^\\p Y \ cos9 BY ) . 

The value of 

w = w{n s E* B E* D «-\p*\\p^\ cos (3* 
m B m D . 

cannot be determined since the angle j3* between the B 
and the D*° in the e + e~ rest frame is unknown. However, 



(3* is bound between a minimum and a maximum value, 
and 

to = [w((3* min ) + w((3* max ))/2 

is a good estimator for w in each event, to and to span a 
range from 1 to 1.5, the distribution of to — to is nearly 
Gaussian with an RMS of 0.026. We use 10 equidistant 
bins of to, their width corresponds to about 2 RMS. 

The fit for the two parameters V = F(l)\V c b\ and 
p 2 is a binned maximum-likelihood fit in three dimen- 
sions. The fit function is the sum of the expected signal 
function S = S(V,p 2 ) and the various expected back- 
ground functions. The signal function in each to bin is 
taken as the product of one-dimensional functions of Am 
and cos 9 BY . These two functions are obtained from fits 
to the reweighted signal MC distributions with V-, p 2 -, 
i?i(l)-, and ^(Independent weights on the generator 
level. The signal fit function also includes the normal- 
ization to the total number of 226 x 10 6 produced BB 
pairs, all decay fractions of sequential decays, the B~ 
lifetime, all MC reconstruction efficiencies, and efficiency 
corrections derived from control data samples and their 
MC expectation. Efficiency corrections for track recon- 
struction and charged particle identification follow those 
of other recent BABAR analyses. For the correction of 
the 7T° reconstruction efficiency we use a control sample 
of r-lepton decays as described below. Small corrections 
are also applied for deviations of the shapes of the Am 
distributions in data and MC because of track resolu- 
tion differences, and for deviations in the shapes of the 
cos 9 BY distributions because of storage-ring energy cal- 
ibration and resolution. 

The background expectation functions are separately 
determined for 23 classes of backgrounds. This large 
number of background functions was necessary in or- 
der to express each function B\ $ as the product of 
Bi^rofAm) and i?2.i,i«(cos 0g Y )- The onc-dimensional fit 
functions are again obtained from fits to MC dis- 

tributions. The fit to the data has 49 free parameters, 
in addition to V and p 2 there are 47 for adjustments of 
background normalizations and shapes, Am shapes, and 
cos9 BY shapes. 

Before fitting the expectation function to the data, it 
is fitted to five different MC subsamples whose size corre- 
sponds to the one of the data sample. All five results for 
V and p 2 agree with the MC input by better than one 
standard deviation. When applied to the data, the fit 
result is V = (36.32 ± 0.60) • 10" 3 and p 2 = 1.146 ±0.055 
with a correlation coefficient g = +0.90. Integrating 
dB/dto over all to leads to B = (5.71 ± 0.08)%. The 
total number of signal events is found to be 23 499 ± 329. 
Though the fit is maximum-likelihood, a control value 
of x 2 can be calculated after the fit as a goodness-of- 
fit measure. We find 4436.3 for 4095 degrees of freedom 
which is, purely statistically, 3.8 a too high. Inspecting 
the distribution of per-bin contributions to x 2 m a U bins 
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TABLE I: Summary of input parameters. 



Input Parameter Value 



Ref. 



B (T(4S) -> B+B- 
B D*° — D°tv°) 
Bm -> K-n+) 
B (tt -> 77) 

T B - 
fll(l) 

^2(1) 



(50.6 ± 0.8)% [12] 

(61.9 ± 2.9)% [12] 

(3.80 ± 0.07)% [12] 

(98.798 ±0.032)% [12] 

(1.638 ± 0.011) ps [12] 

1.417 ±0.075 [6J 

0.836 ± 0.043 [6J 



of ?D, Am, and cos#b Y , we nn d n0 concentrations of high 
values in any area. 

Figure shows the result of the fit together with the 
selected data. The "Signal" part of the fit function 
contains the correctly reconstructed B~ — > D*°e~V e 
decays. The two D** parts contain B — > D**ev de- 
cays with ("Am peaking") and without ("Am flat") a 
correctly reconstructed D*° intermediate state (£)** = 
D\, Dq, D'i, £>2, -D*7r, Dn). Events with a correctly re- 
constructed D*° and a correctly identified electron from 
the same B and from two different B mesons are in the 
"Correlated" and "Uncorrelated" background parts, re- 



spectively, 
and B° - 



"Signal like" arc true decays B 



D 



e v e 



D 



which arc not correctly recon- 



structed. The background from true B — > D°ev decays 
is called u D°ev" . All other background candidates from 
BB events ("Combinatorial D*°") are flat in the Am 
and the cos #g Y distribution since they do neither contain 
a correctly reconstructed D*° nor do they come from a 
charmed semileptonic decay. The last contribution, only 
visible in the highest w bins, comes from cc events in the 
continuum. 

The systematic uncertainties are divided into analysis- 
internal and analysis-external ones, see Table [Til The 
former are specific to our analysis, the latter enter by in- 
put parameters taken from other measurements. Starting 
with the internal ones, the relative uncertainty on the ef- 
ficiency to find a charged particle's track is 0.8%, leading 
to 2.4% and 1.2% for B and V. The dependence of the 
tracking efficiency on the transverse momentum px has 
an uncertainty which could distort the shape of the w 
spectrum. The uncertainties arising from the identifica- 
tion of charged tracks as electrons or as kaons contribute 
to the result as listed under "particle ID efficiency". A 
significant fraction of the total uncertainty of our result 
comes from the precision of the ir° reconstruction effi- 
ciency (e^o). It is determined from e + e~ — ► t + t~ events 
where one of the two r leptons is either reconstructed by 
one track and two clusters (mainly t — > p(-KTt°)v) or it is 
reconstructed by only one track without clusters (mainly 



r — * ttv^vv) llj. The other r, used as a r-pair tag, is 
reconstructed in the channel evv. From the numbers of 
t + t~ events reconstructed in each of the two channels 
we derive an efficiency in data and in MC, giving a cor- 
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FIG. 1: Data distributions (dots 
fit results (stacked histograms) for ( 



with error bars) 
1) Am in the cot 



and 



signal range (-1,±1), (b) cos#b Y i n the Am signal range 
(140,144MeV/c 2 ), and (c) w in both signal ranges. The plot 
below (c) shows the ratio fit/data. The different contributions 
to the fit function are explained in the text. 



7 



TABLE II: Summary of relative systematic uncertainties in 
percent. 





AV/V Ap 2 /p 2 


AB/B 


tracking efficiency (e tr ) 


1.2 




2.4 


Pt dependence of etr 


0.3 


0.5 


0.2 


particle ID efficiency 


0.9 


2.0 


1.6 


extrapolated tv° efficiency (e x o) 


1.8 




3.6 


p n o dependence of e v o 


1.0 


3.5 


0.4 


Am shape of D** background 


0.1 


0.1 


0.2 


shape parameters 


1.0 


2.5 


0.6 


number of BB events 


0.6 




1.1 


off-peak luminosity 


0.1 


0.4 


<0.1 


total internal 


2.8 


4.8 


4.8 


Ri(l) and #2(1) 


0.1 


4.7 


0.3 


f2 {'VIA Q\ . + ^ 
O [1 J — ► -D D f 


0.8 




1.0 


B \d*° -» D%°) 


2.3 




4.7 


B(D° K~n+) 


0.9 




1.8 


B~ life time 


0.3 






D** decay fractions 


0.3 


0.7 


0.3 


number of D*° in cc events 


0.2 


0.7 


<0.1 


total external 


2.6 


4.8 


5.3 


total 


3.9 


6.8 


7.2 


TABLE III: Derivatives of V, 


p 2 , and B. 




V 


P 2 


B 




d/dRi -0.00038 H 


-0.0303 


+0.00382 




d/dRz -0.00158 


-1.22 


+0.00551 





rection to the simulated 7r efficiency. The correction is 
obtained for momenta above 350MeV/c and has a pre- 
cision of 3%. We add to this value 2% in quadrature, 
which is our uncertainty estimate for the extrapolation 
to the lower-momentum range with all ir° mesons from 
D*°ev decays. From fit results for different cuts on p„o 
we estimate the uncertainty in the shape of the w spec- 
trum which gives one of the major contributions to the 
uncertainty of p 2 ("p ff o dependence of e^o"). Corrections 
to the Am shape and to the cos 6*g Y shape are described 
by a paramctrization of the w dependence which also 
contributes to the final uncertainties, see "shape param- 
eters". The determination of the total number of BB 
events in the analyzed data sample has a relative uncer- 
tainty of 1.1%. It contributes only to V and B but not 
to p 2 . The uncertainty on the luminosity of the off-peak 
data sample propagates also to the final result. 

The dominant contribution to the external uncertainty 
on p 2 comes from R\{1) and i?2(l)- We determine the 
derivatives of our fit result with respect to R\ and R%. 
Wc find the values given in Table Hill leading to the uncer- 
tainties listed in Table HT1 The input decay fractions only 
contribute to V and B. An improvement in the precision 
of B (D*° — ► Z?°7r°) would significantly improve our re- 
sults on V and B. The uncertainty on the lifetime t b - 




FIG. 2: One-standard-deviation contours (stat. and sys. 
combined) of V vs. p 2 [l[ from all recent B° — > D* + e~V e 
results together with our preliminary result. The lines cor- 
respond to constant decay fractions B (B~ — > D*°e~77 e ) = 
5.71% ± lcr and ±2cr with a = 0.42%. 

affects only V. Scmilcptonic B decays into higher excited 
charmed mesons, B — > D**ev, contribute to the final un- 
certainties mainly due to their less precisely known decay 
fractions but also due to their description in the fit. The 
uncertainty in the number of correctly reconstructed D*° 
mesons in e + e~ — ► cc events influences B by less than 
0.1%. 

Adding all systematic errors in quadrature leads to the 
last line in Table [TT] and to our preliminary results 

F(l) ■ \Vcb\ = (36.3 ± 0.6 ± 1.4) • 10" 3 , 
p 2 Ai = 1.15 + 0.06 + 0.08 , 
B(B~ -> D*°e-V e ) = (5.71 ± 0.08 ± 0.41)% . 

The correlation coefficients between F(l) ■ \V c b\ and p\ 
arc g = +0.90 for statistics, +0.43 for systematics, and 
+0.52 in total. 

Using F(l) = 0.919 ± 0.033 from lattice QCD we 
obtain \V cb \ = (39.5 + 0.6 + 2.0) • 10" 3 in good agreement 
with the average from the exclusive neutral B decays 
5° -> D*-£+p, (39.2 ±0.7+ 1.4) • 10" 3 Q], and in agree- 
ment with results from the inclusive decays B — > X c £v, 
e. g. (42.0 ± 0.2 ± 0.7) • 10" 3 in Ref. Q- Our result 
for p 2 is in the center of the range (0.5, 1.5) from the 
B" -> D*~e+v experiments Q. 

Figure [2] shows the lcr contour of our result in the 
p 2 , F(l)\ V c b\ plane together with the contours of the 
neutral- B decays. The quasi-diagonal lines in this Figure 
are lines of constant decay fraction B. 

For a comparison of our decay-fraction result with 
the decay fraction of the neutral- B decay mode, we use 
the lifetime ratio t(B+)/t(B°) = 1.076 + 0.008 and 
B(B a D*-£+v) = (5.28 + 0.18)% fl. From this, wc 
expect B(B- -> D*°e~V) = (5.68 ± 0.20)%, again in 
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good agreement with our result. On the other hand, our 
decay-fraction result is about 1.6 a lower than the PDG 
results from CLEO and ARGUS 



121 of the B~ 



average 
0,1- 
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